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Abstract 

The Drude and superfluid weights are investigated in the context of vari- 
ational theory. It is shown that a commonly applied procedure to obtain the 
Drude weight actually corresponds to the superfluid weight. This realization 
changes previous interpretations of conductivity for variational wavefunc- 
tions, one result being that the Gutzwiller wavefunction displays a finite su- 
perfluid weight at finite interaction. To gain insight into DC conductivity an 
expression for the current response in terms of a geometric phase is derived. 
Applying to the Gutzwiller wavefunction it is shown that for any value of 
the interaction the Gutzwiller wavefunction is a perfect conductor. 

1 Introduction 

Variational studies have contributed greatly to our understanding of corre- 
lated systems. In part this is due to their relative simplicity, applicability 
to larger sizes irrespective of the number of dimensions, and their easily 
accessible physical insight. In the case of the Hubbard model [Hill [3] Si fre- 
quently used variational wavefunctions include the Gutzwiller dEl (GWF) 
and Baeriswyl wavefunctions (5j|6l|7l (BWF). The former is based on sup- 
pressing charge fluctuations in the noninteracting solution, the latter on pro- 
jecting with the hopping operator onto a wavefunction in the large interac- 
tion limit. 
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The GWF has been studied by a variety of methods. It can be solved 
exactly in one (U and infinite ifTOl |9j [HI dimensions, and it can be sim- 
ulated in two and three dimensions by variational Monte Carlo lfT2l . The 
one-dimensional exact solution produces a state with a finite discontinuity 
of the momentum density at the Fermi surface. Millis and Coppersmith [ 13 ] 
have investigated the response of the GWF and have concluded that it is 
metallic with a conductivity proportional to the kinetic energy. Insulating 
behavior in projected wavefunctions similar to the one due to Gutzwiller 
can be produced by generalized projection operators Ifl4l [T5l . for example 
non-centro-symmetric or singular projectors. 

Calculating the Drude or the superfluid weights in a variational context 
is a difficult issue. These two quantities can be cast in terms of identical 
expressions (see Eq. (Q])), the second derivative of the ground state energy 
with respect to a Peierls phase [T6l[T9l[T7llT8l . As pointed out by Scalapino, 
White, and Zhang, the two quantities differ in the interpretation of the deriva- 
tive lfT7l[T8l . For the Drude weight the Peierls phase shifts the ground state 
energy adiabatically, remaining always in the same state, for the superfluid 
weight level crossings are also considered. 

In this work general expressions for the Drude and superfluid weights 
are derived in a variational setting. For the Drude weight deriving an easily 
applicable expression is a difficult issue, since the expression derived herein 
depends on the exact eigenvalues of the perturbed Hamiltonian, in practical 
settings often not available. It is then demonstrated that the linear response 
expression for the current can be cast in terms of a geometric phase. The 
tools for calculating this geometric phase (the total position shift operator) 
are also presented. The formalism is then used to interpret the current re- 
sponse of projected wavefunctions. It is demonstrated that the current re- 
sponse in the commonly used Gutzwiller and Baeriswyl projected, as well 
as wavefunctions based on combinations of the two projections, produce a 
current response identical to the wavefunction on which the projections are 
applied (the Fermi sea or the wavefunction in the strongly interacting limit). 

2 Drude and superfluid weights in variational 
theory 

An expression for the frequency (to) and wave vector (q)-dependent con- 
ductivity was derived by Kohn ifToll . The DC conductivity (Drude weight, 
D c ) corresponds to the strength of the ^-function peak of the conductivity in 
the zero frequency limit. The correct expression for D c is obtained by first 
taking the limit (q — > 0) and then the other limit to — > 0. D c is often ex- 
pressed HU [19l in terms of the second derivative of the ground state energy 
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with respect to a phase associated with the perturbing field as 

~d 2 E (®Y 



D '-L 



<9<D 2 



(1) 



O=0 



Here Eq{<&) denotes the perturbed ground state energy, CD denotes the Peierls 
phase. 

Scalapino, White, and Zhang (SWZ) IfTTl [T8 ll have investigated the dis- 
tinction between the Drude and superfluid weights. In particular they studied 
the importance of the order of different limits (a> — > 0, q — > 0) for the con- 
ductivity. In a variational context implementation of the frequency limit is 
not straightforward, since, strictly speaking there is no frequency to speak 
of. However, SWZ have also pointed out that the derivative with respect to 
the phase O in Eq. £[]) is ambiguous. They showed that if the derivative is 
defined via adiabatically shifting the state which is the ground state at zero 
field, then the Drude weight results. In the presence of level crossings the 
adiabatically shifted state may be an excited state for a finite value of the 
perturbation. The superfluid weight is obtained if the derivative corresponds 
to the "envelope function", i.e. the ground state of the perturbed sytem is 
taken to define the derivative. The distinction between these two derivatives 
can be implemented by embedding the periodic system under study in a 
larger periodic system, and defining the perturbation in terms of the periodic 
boundary conditions of this larger system. In cases in which level crossings 
are close to (3> = 0) conductors, superconductors, and insulators can be dis- 
tinguished IfTTl [T8lf20ll . In general, the position of level crossings depends 
on dimensionality IfTTl fT8ll . 

The finite temperature extension of D c has been given by Zotos, Castella, 
and Prelovsek[21 ] (ZCP). This generalization can be summarized as 



d 2 E n {<S>) 



<9<D 2 



(2) 



o=o 



Note in this expression the Boltzmann weight factors remain unchanged as 
the perturbation <I> is turned on. Eq. (© has been applied If22l to calculate the 
DC conductivity in strongly correlated systems. Taking the zero temperature 
limit reproduces Kohn's expression for D c . To define a finite temperature 
analog of D s one lets the Boltzmann weight factors depend on the perturbing 
field CD as 



DAT) = I 



d 1 y exp(-/3£ n (Q)) 
Q 



(3) 



o=o 



Indeed the ground state superfluid weight is reproduced in the zero temper- 
ature limit. Eqs. © and ([3]> follow from the assumption that the distinction 
between the Drude and superfluid weights is due to the different types of 
derivatives as discussed by SWZ. 
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Similarly, in deriving expressions for D c and D s in a variational setting 
our starting assumption will also be that the distinction between the two 
quantities is due to the effects of level crossings. Suppose is a varia- 

tional wavefunction, where y denotes a set of variational parameters, which 
we wish to use to optimize some Hamiltonian H with eigenbasis 

H\W n ) = EM. (4) 

The estimate for the ground state energy may be written in terms of a density 
matrix as 

<*( 7 )|d|*(y)> - 2 <t ^ )|VF " >£ "< VF " l ^ )> = Z PnE >" (5) 

n n 

the probabilities can be written as 

P n = l<*(y)|T„>| 2 . (6) 

Comparing with Eq. © it is obvious that a consistent formalism requires 
that the variational Drude weight be defined as 

n 

with P„ independent of the perturbation cD. It follows that the variational 
parameter y is independent of the perturbation The variational analog of 
D s (based on Eq. (O) corresponds to 

n 

where the probabilities P n (<£>) depend on <I> and the variational parameters y 
in this case also depend on O. 

A central difficulty in calculating D c in a variational theory is the fact 
that it depends on the exact eigenvalues of the perturbed Hamiltonian (see 
Eq. ©), however variational theories are usually applied in cases where the 
exact solution is not easily accessible. While the Drude weight remains a 
difficult problem in general, it is shown below that the current in general can 
be cast in terms of a geometric phase, and evaluated even in a variational 
context. 



3 Current in terms of a geometric phase 

In this section we consider the adiabatic current response of a system in gen- 
eral, not only in a variational context. After showing that the persistent cur- 
rent can be expressed as a geometric phase [23] EH, we explicitly construct 
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the mathematical tools to calculate it, and use the results in the next section 
to interpret the GWF. Since the current can be cast in terms of observables, 
it follows that the calculation of the Drude weight is also accessible, being 
the first derivative of the current as a function of the Peierls phase. 

Consider a system periodic in L, and experiencing a perturbation in the 
form of a Peierls phase cp. Its Hamiltonian can be written as 

" (A + <D) 2 



6(0) = J] KPi 2m + V(x u x N ). (9) 
The following identity also holds 



e^m = V . (io) 

m 

r=l 

The ground state energy can be written as 

E(0) = <¥(®)|^(0)|¥(®)>. (11) 
The average current for such a system can be expressed as |[T6l 

/(©) - <9®£(<D) = <T(0)|d<D^(0)|Y(0)>. (12) 
Substituting for the partial derivative of the Hamiltonian we obtain 

Jm = — + ywm-m<i>))- (i3) 

i=i 



In the position representation the current can be written 

Jm = V W) | |vp ( 0)). (14) 

m m j—* oxi 

Next we rewrite the wavefunction in terms of a shift of the total position and 
define a wavefunction 

<jci, ...jc*PP(G; X)> - T(jci + X, x N + X; O). (15) 

The action of the total momentum can then be cast in terms of the derivative 
with respect to X as 

V — ¥(jci + X, jcat + X; <D) - d^Oq + X, x w + X; O). (16) 

The effect of X on the particle positions is similar to the effect of the Peierls 
phase on the momenta. Like the Peierls phase it is an external parameter, so 
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one can perform adiabatic cycles as a function of it. Expressing the current 
in terms of averaging in X over the unit cell £ J () dX... leads to 

N® i r L 

7(0) = dX<T(0;X)|5 x |T(0;X)>. (17) 



The second term in Eq. (fPTT ) is a geometric phase [23]. Since it results 
from the periodicity of the parameter X it is similar to the geometric phase 
derived by Zak |[25l . It is also similar to the geometric phase expression 
appearing in the modern theory of polarization l26l . with the variable X 
playing the role of the crystal momentum in this case. Thus the current due 
to a perturbation can be expressed in terms of a constant proportional to the 
number of particles, and a geometric phase term. Below an interpretation of 
the phase term is given. It is interesting to note that a finite persistent current 
is in principle possible for an unperturbed system (the case O = 0). 

The next question to address is the actual calculation of this quantity. We 
can construct a scheme which is in the spirit of the total position operator 
proposed by Resta ETl 1281 to calculate the polarization. We consider the 
case O = (and suppress the notation), without loss of generality. We first 
rewrite the Berry phase appearing in the expression for the current in terms 
of its discretized analog as [29 ] 

. M-1 

7(0)= lim — Imln n< , J , (^AX)| l P((5+ 1)AX)). (18) 
ax^o mL 1 1 
s=0 

The continuous expression can be recovered by Taylor expanding the wave- 
function + l)AX)) around sAX and taking the limit as AX — > 0. Indeed 

j M-1 

7(0) = lim — Im V lnOF(jAX)| [|¥(s)AX)> + d x msAX))AX] = ln[l+( l P(sAX)|5xl x I 7 (iAX)>AX]. 
ax^o mL t—i 

.v=0 

(19) 

When the limit AX — > is taken the natural logarithm can be expanded and 
we obtain 



7(0) = — Im 
mL 



J dxmx)\d x mx)) = ~ f dx(y(x)\d x mx)). 

(20) 

The shift in the total position of the system by a value of AX can be accom- 
plished using the total position shift operator U(AX). The explicit form of 
this operator will be derived below, for now we assume its existence. We 
define it as 

U(AX)\ x l'(X)) = \Y(X + AX)), (21) 
Using Eq. (1211) we can express the product in Eq. (TT8l) as 

M-1 

Y\ (^AX)!^ + 1)AX)> = ( V F(0)|?7(AX)|»P(0)) M (22) 

s=0 
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Substituting into Eq. (fT8l > the expression for the current becomes 

7(0)- lim -— Imln^OWCAX^CO)). (23) 
ax^o m AX 

The total position shift operator can be constructed using real space per- 
mutation operators. This derivation has been given elsewhere [30], here we 
emphasize the main results. In second quantized notation the permutation 
operator between two positions can be written as 

P ij = l-(c] -c))(ci-Cj). (24) 

This operator has the properties 

Pifi) ~ CiPij'Pifii ~ c jPij>Pij c ) = C i^ij^ij c i = C )^'j- @5) 

Assuming a grid with spacing AX, using Pjj we can construct an operator 
which shifts all the positions on the grid in a periodic system. The operator 

U(AX) = P 12 P23...Pl-il, (26) 

where it is assumed that the indices refer to particular grid points, has the 
property that 

A fAy , / c^U(AX), i = 2,...,L 

U(AX)c, = l cl0(AXI . = l (27) 



It also holds that 



U(AX)c k = e iAXk c k U(AX), (28) 



where c k denotes the annihilation operator in reciprocal space. Eq. (1281) can 
be demonstrated by Fourier transforming c k and applying (1271) . Taking the 
Fermi sea 



as an example one can show that 



U(AX)\FS ) = e iAXK \FS ), (30) 



with K = h 



As an example we consider again the non-interacting Fermi sea given by 

ire> = g£.4 r IO>. (3D 

where the ^-vectors are spread symmetrically around zero. Applying a per- 
turbation <I> shifts all /c-vectors by O. The resulting current is 

2N 

7((D) = _o, (32) 
m 
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corresponding to a Drude weight of D c - 2N/m. It is interesting to see that 
the current is proportional to twice the number of particles. In a Fermi sea 
conduction can occur due to particles as well as holes, of which at half-filling 
there are an equal number. For systems with bound particles and holes, /(cD) 
is reduced, as bound excitons do not participate in conduction and reduce 
the effective number of charge carriers. Thus the geometric phase in Eq. 
(fTTl) accounts for exciton binding. When all the particles are bound to holes 
then the constant term in Eq. (TTTT ) is cancelled by the phase term leading to 
/(<!>) = 0. An example of bound particles and holes in the same band leading 
to insulating behaviour is the Baeriswyl variational wavefunction O. 



4 Current response of projected wavef unc- 
tions 

In this section we provide the response theory of some commonly used pro- 
jected wavefunctions HJEHUIHl. We emphasize that it is the current response 
we calculate, not the Drude weight, which is the first derivative of the current 
response with respect to the perturbing phase. 

The Gutzwiller wavefunction HH3 (GWF) was proposed as a variational 
wavefunction for the Hubbard model, and it has the form 

[¥ G {y)) = e^ A \FS), (33) 

where D - Yui n i\ n ii- Without loss of generality we consider the one- 
dimensional case. 

Before developing the current response theory of the GWF, we present 
the calculation of a quantity which expresses the extent of localization. Lo- 
calization has been suggested long ago as a general criterion of metallic- 
ity l[T6l . and the relation of the spread to the DC conductivity has been 
shown in a number of places ETl 1281 l20ll . In particular we calculate the 
normalized spread defined as 

(34) 

Due to the ill-defined nature of the position operator in periodic systems we 
choose a sawtooth representation which can be written as 



X 



L/2-1 ( l p„ 



■2 



m=-L/2 

m±0 



2 + exp ( i2f ) . 



(35) 



where W denotes the total momentum shift operator, which has the property 
that 

WmK)) = \¥{K + (2tt)/L)>. (36) 
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L 


Fermi sea 


y = 1.0 


y = 2.0 


12 


0.08275 


0.079(1) 


0.0412(9) 


24 


0.08312 


0.0830(6) 


0.0682(6) 


36 


0.08327 


0.0831(5) 


0.0797(5) 


48 


0.08330 


0.0833(4) 


0.0824(4) 


60 


0.08331 


0.0829(3) 


0.0830(3) 


oo 


1/12 







Table 1: Spread in the total position divided by the square of the system size 
for the Fermi sea and the Gutzwiller wavefunction. Two different values of the 
variational parameter, y = 1 .0 and y = 2.0 are shown. As the system size increases 
the value 1/12 is approached by all three systems. The approach to the limiting 
value slows down as correlation effects are introduced, it is slowest for y = 2.0, 
the "most projected" of the three examples. 



The construction f3Tj of this operator is analogous to the total position shift 
operator used to define the persistent current in section [5] For a state |cD) 
diagonal in the position representation one can write 

ty\q>) = J%2'*'\Q) (37) 

where xt denotes the position of particle i. Using the sawtooth representation 
one can show that for the Fermi sea that the spread in position 

(X 2 )-(Xf 1^ 1 1 
Ll "^2(l-cos(^))-12- 

To show this one needs to substitute Eq. (T35T ) into Eq. (1341 ). and then use the 
fact that 

W'"\FS) = i ° m = 2 — L - 1 (39) 
1 ; \ \FS) m = 0. K ' 

Our results are shown in Table [TJ The GWF results for two different 
values of the variational parameter were calculated for a one-dimensional 
system based on the variational Monte Carlo method of Yokoyama and 
Shiba 02). The fact that the normalized spread approaches a constant for 
large L (system size) indicates that the system is delocalized, hence metallic. 
What is surprising in these results, however, is that for large L the spread of 
all three examples converges to the same value. The projecting out of double 
occupations in the GWF seems to have no effect on the spread for large L, 
and is identical to the result for the Fermi sea. The GWF though is thought 
to be a representative of "bad metals", metals whose conductivity is reduced 
due to strong correlations. |[32l[33ll 
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It turns out that these results are actually consistent with what one ob- 
tains for the current response. We consider the phase term under a perturba- 
tion in the form in Eq. (1231) . Consider first the action of the operator U(AX) 
on the GWF. 

#(AX)PF(y)) = U{AX)e' y ^' n '^\FS). (40) 

The operator U (AX) shifts the positions of all particles by one lattice spac- 
ing. Such a shift will not change the total number of double occupations, 
hence the Gutzwiller projector and the total position shift operator commute. 
We can write 

U(AX)\¥(y)) = e- yZ ' n '^U(AX)\FS) = e'^'^y)), (41) 

where kj denote the sum over the momenta of the Fermi sea. One obtains 
exactly the same result in the absence of the Gutzwiller projector. When 
substituting back into Eq. (1231) we find that the current response of the GWF 
is exactly that of the Fermi sea, and this result is independent of correlations 
(whose strength increases monotonically with the variable y). The above 
derivation can be extended to projections based on Jastrow-type correlations 
and the conclusion is valid as long as the projections are centro-symmetric 
(considered in Ref. ifTBlQ . It has been shown [ 14] that non-centrosymmetric 
correlators can produce an insulating state. The current response in this case 
will also not follow the derivation above, since a shift in all the particles can 
change the contribution to the projector. Another exception is the case when 
y — > oo, i.e. the singular case, which in general is also known to allow for 
insulating behavior lfT3l . 

For the GWF one can obtain further insight into the current response by 
writing it in the position representation as 

|*F G > = ^ e~ yD(R) Det(K; R)|R>. (42) 

R 

In Eq. (1421 R indicates the configurations of particles (both up-spin and 
down-spin), D(R) indicates the number of double occupations for a par- 
ticular configuration of particles, Det(K; R) denotes the product of Slater 
determinants for up-spin and down-spin electrons, and |R> denotes a posi- 
tion space eigenstate. From Eq. (1421 . we see that the projection changes the 
relative weight of different configurations but leaves their phases intact. 11331 
The fact that the current, a quantity related to the phase of the wavefunction, 
is unaltered by the Gutzwiller projection coincides with the result above, 
that the persistent current for a Gutzwiller wavefunction is determined ex- 
clusively by the Fermi sea. In fact Millis and Coppersmith [13 ] suggest a 
scheme in which a projector operator of the form e' s , with S = jj(Hf -H~), 
(Hf(H~) raises(lowers) the number of double occupations) acts on the Fermi 
sea to produce an insulating wavefunction. Clearly this scheme would alter 
the phases of the Fermi sea. 
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The above reasoning can be extended to other commonly used projected 
variational wavefunctions. The Baeriswyl-Gutzwiller wavefunction can be 
written 

\y B G(a,7)) = e- af e-? 6 \FS). (43) 

In this case T denotes the kinetic energy, and a denotes the variational pa- 
rameter. Since the total position shift operator U(AX) is diagonal in momen- 
tum space, it trivially commutes with the projector e~ oT . We conclude that 
the current response of the Baeriswyl-Gutzwiller projected wavefunction is 
identical to that of the Fermi sea. The other two commonly used variational 
wavefunctions are the Baeriswyl and Gutzwiller-Baeriswyl projected wave- 
functions. Their form is 

PF B (ar,y)> = e-^PFoo), (44) 

\y G B(<*,y)) = e-y 6 e- at \y co ). (45) 

In Eqs. (1441 and (|43T > denotes the wavefunction in the limit of infinite 
interaction. This function is in general not known. Again one can exploit the 
fact that the total position shift commutes with the projector operators and 
conclude that the current response in these both cases will depend on l^oo) 
exclusively. While this functions is not known, in general, in the half-filled 
case one can assume that its current response is zero. 

5 Conclusion 

The current response was investigated in the context of variational theory. 
The Drude and superfluid weights have seemingly identical expressions (sec- 
ond derivative of the ground state energy with respect to the Peierls phase), 
however, as was pointed by Scalapino, White, and Zhang, the meaning of 
the derivative differs between the two, one being the adiabiatic the other 
the "envelope" derivative. Assuming their interpretation of the derivative 
we derived the appropriate expressions for the Drude and superfluid weights 
appropriate in variational theory. A key difficulty with the former is the 
appearance of the exact eigenstates of the perturbed Hamiltonian, in gen- 
eral not available in practical situations where variational theory is used. As 
a partial remedy the persistent current was shown to consist of a constant 
term, proportional to the perturbation and the number of charged particles, 
and a geometric phase term. This expression can be used in practical set- 
tings to obtain the Drude weight by numerically taking the first derivative of 
the current with respect to the phase. The current response of several com- 
monly used variational wavefunctions was also analyzed, and shown that 
variational wavefunctions which use a Baeriswyl or Gutzwiller projection 
will have a current response determined by the wavefunction on which the 
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projectors are applied (Fermi sea or the solution in the strongly interacting 
limit). 
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